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ABSTRACT 



Using the correlation between the equivalent width of 
the diffuse interstellar band (DIB ) at 8620 A and the in- 
terstellar reddening reported by iMunaril (120001) GAIA 
could directly trace the interstellar extinction through- 
out the Galaxy. We checked for the magnitude and dis- 
tance limitations of this method in 42 Galactic directions 
by simulating the RVS data on the stellar sample pro- 
vided b y the model of stel lar population synthesis of the 
Galaxy (Robin et al. 2003). The simulation indicates that 
the imprint of the 8620 A DIB will be detected in the RVS 
spectra of stars with magnitudes up to V ~ 16 with suf- 
ficient accuracy to trace not only the distribution of the 
interstellar medium but also the radial component of its 
kinematic motion, i.e. the radial velocity of the mass cen- 
ter of the dust cloud in the line of sight. 
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1. INTRODUCTION 



The purpose of the here presented simulation was to 
check for the limitations and to estimate the distance 
range of this interstellar mapping method. 



2. SIMULATION 



The Galactic stellar population was simulated with 
the Besancon model of stellar population synthe- 
sis (iRobinet al J 12003). The interstellar extinction was 
taken into account with the use of t he three-d imensional 
model of the dust distribution by iDrimmel & S pergel 
(1200 ll) and lDrimmel et al.1(l2003l) . Simulations were per- 
formed in 42 Galactic directions (/ = 10°, I = 20°, 
I = 30°, I = 60°, I = 90°, I = 135° and I = 180°; 
b = 0°, b = 2°, b = 5°, b = 10°, b = 30° and b = 60°) 
with the solid angles varying between 0.2 and 15 square 
degrees for different Galactic directions in order to keep 
the number of the simulated stars approximately constant, 
i.e. between 5000 and 15000. Simulated Galactic direc- 
tions are shown in Figure [0 Two exemplary stellar sam- 
ples are presented in Figure |2] 



A surprisingly good correlation between the equivalent 
width (EW) of the diffuse interstellar band (DIB) at 
8620 A and the inte rstellar reddening (Eb-v) nas been 
reported by Munari( 2000). The relation was calibrated as 
E B -v = 2.69 x EW{k) on the base of 37 stars widely 
distributed in Galactic coordinates and distances. How- 
ever, the exact slope of the relation is expected to depend 
on the properties of the interstellar medium and to vary 
with direction and distance. The most general relation 
would thus be 

Eb-v = a(l,b,D) x EW{k) . 
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GAIA will determine Eb-v from the photometric mea- 
surements. If the 8620 A DIB will be resolved from RVS 
spectra, the fine-grid calibration of a(l, b, D) will be per- 
formed and thus the distribution and the chemical com- 
position of the interstellar matter will be obtained. More- 
over, the exact wavelength position of the 8620 A DIB 
will reveal the information on the kinematics of the inter- 
stellar matter. 



Figure 1. The Galactic stellar population was simulated 
in 42 directions. The solid angle of each simulation was 
chosen in such a way that the number of simulated stars 
was kept approximately constant, i.e. ~ 5000 — 15000. 
The sizes of the solid angles that varied between 0.2 and 
15 square degrees are presented in a relative scale with 
the sizes of the black dots. 
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Figure 2. Two exemplary simulated stellar samples in the directions I — 10°, 6 = 2° (left) and I = 10°, b = 30° ( right). 
Above: Absolute V magnitudes (black) and apparent I magnitud e s (gra y) of the simulated stars as a function of distance. 
Below: Interstellar extinction according to yDzimmel & SpereeA ft20()l\) as a function of distance in these two Galactic 
directions. 



The simulation of the DIB detection accuracy was 
performed as follows. For any star given by the 
Besancon mod el the clo sest computed synthetic spec- 
trum (Zwitt er et alJ20 04) from the multi -parameter spec- 
tral grid was taken. The spectrum was then shifted for 
the star's radial velocity and the 8620 A DIB Gaussian 
line, Doppler shifted for between and the shift of the 
star was added. The width of the Gaussian line was taken 
as in lMunaril lEOQO). i.e. FWHM was taken to be 5.5 A 
and the amplitude was determined from EW . Finally, the 
noise appropriat e to the star's magnitude and the RVS in- 
strument (Mun ari et alJl20 03) was added to the synthetic 
spectrum. In addition to the simulated star spectrum the 
synthetic spectrum template was generated, which dif- 
fered from the simulated spectrum for one step in each of 
the parameters (temperature, metallicity. . . ) o f the spec- 
tral gr id. It also differed in the radial velocity (iKatz et al.l 
2004) and possessed no DIB line. The generated template 
was subtracted from the RVS simulated spectrum and the 
8620 A DIB line was resolved from the residuum by fit- 
ting the Gaussian line with three parameters, i.e. EW, a 
and Ao on it. 



3. RESULTS 



ing out from the Galactic plane the 8620 A DIB line can 
be detected in stellar spectra up to distances 10-15 kpc. 
Finally, at very high latitudes the imprint of small inter- 
stellar extinction becomes weak, which again reduces the 
detection of the DIB line. Tracing of the 8620 A DIB line 
varies also with the Galactic longitude. The Galactic map 
with the achieved distance ranges for the interstellar mat- 
ter mapping via the 8620 A DIB line is shown in Figure|3] 
Results of the 8620 A DIB line recovery from the simu- 
lated RVS spectra in six Galactic directions are shown in 
Figures @]and|5] 
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Interstellar extinction is responsible for both, the imprint 
of the 8620 A DIB line in the stellar spectra and the in- 
crease of stellar magnitudes that limits the detection of 
the DIB line. It is not surprising, that the simulation re- 
sults indicate the following. High interstellar extinction 
is a strong limiting factor in the Galactic plane where 
the achievable range of the interstellar matter mapping 
via the 8620 A DIB line is smaller than 5 kpc. When go- 



Figure 3. The Galactic map with achieved distance 
ranges for the interstellar matter mapping via the 8620 A 
DIB line. This method of the interstellar mapping is 
the most efficient at intermediate latitudes, where the in- 
terstellar extinction is already moderate but still high 
enough to leave the significant imprint in the stellar spec- 
tra. 
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Figure 4. Recovery of the 8620 A DIB line from the simulated RVS spectra: Number of simulated stars as a function 
of distance in six Galactic directions (I = 10°; b = 0°, b = 2°, b = 5°, b = 10°, b = 30° and b = 60° ). Black 
color represents all simulated stars, dark gray those for which the DIB line was recovered with an accuracy of at least 
30 km s _1 and light gray those for which the recovery of the DIB line was more accurate than 1 km s" 1 . 
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Figure 5. Recovery of the 8620 A DIB line from the simulated RVS spectra: Number of simulated stars as a function of 
in accuracy in eq uivalent width, i.e. the difference between the simulated and recovered EW of the DIB line. According 
toM unar EW is proportional to Eb-v, where EW = 0.5 A corresponds to Eb-v — 1-35. 



To conclude, simulated RVS data indicate that the 8620 A 
DIB line will be detected in the RVS spectra of relatively 
faint stars with magnitudes up to V ~ 16 with sufficient 
accuracy to trace not only the distribution of the inter- 
stellar medium but also the radial component of its kine- 
matic motion. However, it is necessary to point out that 
the obtained results of the simulation should be further 
constrained by more careful treatment of the S/N issue. 
Other realistic problems such as crowding in the stellar 
field should also be taken into account. 
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